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The experimental valence band photoemission spectrum of semiconductors exhibits multiple satel- 
lites that cannot be described by the GW approximation for the self-energy in the framework of 
many-body perturbation theory. Taking silicon as a prototypical example, we compare experimental 
high energy photoemission spectra with GW calculations and analyze the origin of the GW failure. 
We then propose an approximation to the functional differential equation that determines the exact 
one-body Green's function, whose solution has an exponential form. This yields a calculated spec- 
trum, including cross sections, secondary electrons, and an estimate for extrinsic and interference 
effects, in excellent agreement with experiment. Our result can be recast as a dynamical vertex 
correction beyond GW, giving hints for further developments. 



Photoemission is a prominent tool to access informa- 
tion about electronic structure and excitations in mate- 
rials. Modern synchrotron sources can provide detailed 
insight, thanks to their high intensity and broad pho- 
ton energy range. But the interpretation of the experi- 
mental data is far from obvious, and theory is an essen- 
tial complementary tool. However, ah initio calculations 
typically focus on bulk bandstructure [l|, |2[; thus sur- 
face effects are ignored, and satellites are not included. 
The latter are a pure many-body effect due to coupling 
to excitations of the material. Such many-body effects 
are contained in approaches developed for correlated ma- 
terials [3|, |j| however, these are usually based on mod- 
els with short-range interactions, whereas satellites such 
as plasmons involve long-range effects. Plasmon satel- 
lites have been extensively studied in core-level experi- 
ments [5|- There they can be described by a theoreti- 
cal model where a single dispersionless fermion couples 
to bosons. The resulting exact Green's function has an 
exponential form given by the so-called cumulant expan- 
sion (CE). A Taylor expansion of the exponential leads 
to a well defined quasi-particle (QP) peak followed by 
a decaying series of plasmon satellites at energy differ- 
ences given by the plasmon energy, consistent with ex- 
perimental observations |6l4lC|. In the valence region, 
plasmon satellites are much less studied, though ah ini- 
tio approaches can provide a good starting point. At high 
photoelectron energies the photoemission spectrum is ap- 
proximately proportional to the intrinsic spectral func- 
tion A{uj) = — (l/7r)Im5(cij), where Q is the one-particle 
Green's function. The latter is typically calculated using 
the widely used GW approximation (GWA) [tI. ITTI. ITi| . 
In principle, the GWA contains correlations effects be- 



yond the quasiparticle approximation. However, these 
additional features are rarely calculated due to computa- 
tional complexity and, more importantly, the serious dis- 
crepancies between GWA and experiment (see e.g. |l3l - 
llGJI). The CE has also been used for the homogeneous 
electron gas [17| and simple metals IJ, ll5| , yielding an 
improved description of satellites over GW. Silicon 16 1 
and graphite [l8| were also studied, but no plasmon satel- 
lite series were observed. However, these results are not 
conclusive due to difficulties of interpreting the experi- 
mental data. This leaves a series of important questions: 
(i) do materials generally exhibit intrinsic satellites in 
the valence band region following a cumulant like distri- 
bution, or are the extrinsic plasmon peaks, due to losses 
incurred by the escaping photoelectron, dominant? (ii) 
if such series are seen, how bad are ab initio GW calcu- 
lations, what is the reason for their failure, and (iii) how 
can they be improved? Answering these questions would 
be a crucial step towards a better understanding of cor- 
relation effects in electronic excitations and a predictive 
ah initio approach to photoemission. 

In this work wc focus on plasmon satellites using sili- 
con as a prototypical example. We have obtained valence 
band photoemission data at high photon energy (XPS) 
that constitute a reliable and well resolved benchmark. 
Analysis of the data allows us to elucidate the failure 
of GW in describing the satellites. Then, starting from 
the fundamental equations of many-body perturbation 
theory (MBPT), we show how the failure can be over- 
come by using a decoupling approximation that leads to 
an exponential representation of the one-particle Green's 
function. Together with an estimate for extrinsic and in- 
terference effects, we obtain results for the quasiparticle 



peaks and satellites in excellent agreement with experi- 
ment. Our theoretical results can be expressed in terms 
of a dynamical vertex correction, a powerful basis for fur- 
ther modelling. 

Angular resolved valence photoemission (ARPES) 
measurements were performed at the UHV photoemis- 
sion experimental station of the TEMPO beamline [19[ 
at the SOLEIL synchrotron radiation source. Linearly 
polarized photons from the Apple II type Insertion De- 
vice (HU44) were selected in energy using a high resolu- 
tion plane grating monochromator with a resolving power 
E/ /S.E = 5000. The end-station chamber (base pressure 
10-1° mbar) is equipped with a modified SCIENTA-200 
electron analyzer with a delay-line 2D detector which op- 
timizes the detection linearity and signal/ background 
ratio '20[. The overall energy resolution was better than 
200 meV. The photon beam impinges on the sample at 
an angle of 43°, and photoelectrons were detected around 
the sample surface normal with an angular acceptance of 
±6°. A n-type {No ^2x 10'^^ P atoms/cm^) Si(OOl) 
wafer was cleaned from the native oxide by flash anneal- 
ing at 1100° C after prolonged degassing at 600° C in 
ultra-high vacuum. The silicon surface was annealed at 
300° C to prevent surface etching, and hydrogenated in 
a partial pressure of activated hydrogen about 2 x 10"^ 
mbar for 20 min. The ARPES was measured along the 
S direction. At 800 eV kinetic energy the Si Brillouin 
zone is observed with an emission angle slightly smaller 
than 5°. The measured photoemission map was inte- 
grated over the spectral intensity originated by two Bril- 
louin zones. The Fermi level was obtained by measuring 
a clean Au(lll) surface. The experimental data (crosses) 
are summarized in Fig. [1] One can distinguish the quasi- 
particle peaks between the Fermi level at zero and the 
bottom valence at -12 eV, followed by two prominent 
satellite structures, each at a mutual distance of about 
17 eV, as well as a more weakly visible third satellite 
between -52 and -60 eV. These structures are obviously 
related to the 17 eV silicon bulk plasmon 2l|. 

The exact one-electron Green's function Q is described 
by an equation of motion with the form of a functional 
differential equation 22], 

G = Go + GoVhG + GovG + iGoVc-7-- (1) 

Here Go is the non-interacting Green's function, </? is a 
fictitious external perturbation that is set to zero at the 
end of the derivation, Vc is the bare Coulomb interaction, 
and all quantities are understood to be matrices in space, 
spin, and time. The Hartree potential Vh gives rise to 
screening to all orders. Linearizing Vh with respect to ip 
yields [23] 

G{tit2) = G"H{tit2) + GUtihMh)G{ht2) 
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FIG. 1. (Color online) Experimental XPS spectrum of Si at 
800 eV photon energy (blue crosses), compared to the the- 
oretical intrinsic A{u}) calculated from GoWo (red dashed), 
and from Eq. (J4]) (green dot-dashed). On top of the latter 
the black solid line also includes extrinsic and interference ef- 
fects. All spectra contain photoabsorption cross sections, a 
calculated secondary electron background and 0.4 eV Gaus- 
sian broadening to account for finite fc-point sampling and 
experimental resolution. The Fermi energy is set to eV. 



where (p is equal to ip screened by the inverse dielec- 
tric function, W is the screened Coulomb interaction, 
and G^ is the Green's function containing the Hartree 
potential at vanishing (p] only time arguments are dis- 
played explicitly and repeated indices are integrated. 
This linearization preserves the main effects of W and 
hence of plasmons. With the additional approximation 
g-u) — G{tzti)G{tit2) one obtains the Dyson equa- 
tion G = Gh + Gh^G in the GWA for the self-energy E. 
However this approximation can be problematic. For the 
following analysis we use the standard GqWq approach, 
where Go is taken from an LDA calculation and Wo is 
the screened interaction in the Random Phase Approx- 
imation. Fig. [2] shows the GqWo spectral function [2J| 
A{lj) = i |ImS(u;)|/[[c^-£ff-ReE(aj)]' + [ImE(c^)]'] 
of Si at the F point, for top valence (solid line) and 
bottom valence (dashed), respectively. The top valence 
shows a sharp quasiparticle peak followed by a broad, 
weak satellite structure at about -21 eV. This peak stems 
from the prominent peak in Im E (full circles) at about 
— 18 eV, itself due to the plasmon peak in ImW. It 
is a typical plasmon sateUite, though (cf. [7|), the QP- 
satellite spacing is slightly overestimated because the 
term cj — eh — ReE (full squares) in the denominator 
of the expression for A{uj) is not constant. However the 
GWA has a more severe problem: for the bottom valence, 
the satellite structure at about —36 eV is much too far 
from the QP peak at about —12 cV, and much too sharp. 
This satellite does not correspond to a plasmon peak in 
ImE (empty circles), but to a zero in w — sh — ReE 
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FIG. 2. (Color online) GoWo spectral function of bulk silicon 
for the top and bottom valence bands at the F point (black 
solid and blue dashed, respectively). The corresponding imag- 
inary parts of the self-energy (red empty-circles-dashed-line 
and green full-circles-solid-line) and Lj — eH— ReE (red empty- 
squares-dashed-line and green full-squares-solid-line) are also 
shown. The Fermi energy is set to eV. 



(empty squares) in the denominator of A{uj), as for a 
QP peak. It has been interpreted in the HEG as a plas- 
maron, a coupled hole-plasmon mo de 1251 . but as noted 
below it is an artifact of the GWA [261 127| . Fig. [1] com- 
pares the total GW spectral function (dashed red line) 
summed over all valence bands and fc-points, with our 
XPS data. The effects of cross-sections are included by 
projecting on angular momenta in atomic spheres using 
the atomic data of Ref. 



28|. The secondary electrons 



background at energy uj was determined by integrating 
the calculated intrinsic spectral intensity between u and 
the Fermi level, similar to [29j. A constant scaling factor 
was set such that the measured photocmission intensity 
at the highest binding energy (60 eV), where primary 
electrons intensity is absent, is reproduced. As expected, 
the dominant QP spectrum is well described by GW, but 
the satellite is dominated by the plasmaron around —36 
eV, in complete disagreement with experiment. The ex- 
perimental plasmon satellite at about —25 eV appears 
only as a weak shoulder in the GWA. Thus the plas- 



26|, l27| 



maron peak is responsible for the GWA failure 
in silicon. 

To go beyond the GW self-energy requires vertex cor- 
rections. However, adiabatic vertex corrections (see e.g. 
|30| ) only lead to renormalization of energies and do not 
create new structures. Thus alternatively, we concen- 
trate here on dynamical effects, and we choose to ap- 
proximate directly Eq. ([2]), without passing through a 
self-energy. We decouple Eq. (^ approximately by sup- 
posing that Q and Qh are diagonal in the same single 
particle basis. Eq. ([2]) is then applicable separately for 
every single matrix element of Q and each state couples 



independently to the neutral excitations of the system 
through W [311]. The latter can now be understood as 
the screened intra-orbital Coulomb matrix element for 
the chosen state. Such a decoupling approximation can 
be optimized [26|, |27[ by adding and subtracting a self- 
energy correction, hence by using a QP Green's func- 
tion Qa obtained from a good QP self-energy instead 
of Qh- Since the GWA is currently the state-of-the 
art for QP properties, we suppose that for every decou- 
pled state k, G^i''') = i9{—T)e~'^^'''^ is determined from 
SG^(efc), where e^ = 6° + ^^^{sk) is the (complex) 
GW quasiparticle energy and r = ti — ^2- Now Eq. ([2]) 
can be solved exactly for each state. Briefly the main 
steps are: (i) solve the non-interacting (W = 0) ver- 
sion of (I2|), which leads to an explicit solution Q'^; (ii) 
iterate the result g ^ g^ ~ g^Ag + ig^W^ start- 
ing from tj(°) = g'^. Here A compensates for the self- 
energy insertion used for the optimized decoupling; (iii) 



use the exact relation 



<5¥'(*4) 



iQ'i{Ut2)6{t2 - t4)6l(i4 - ti) to derive 



Q{tit2) - ^AMe'^^e 



Ar i/,'f df[(p(t')-Jl^ dt"W(t't")] 



(3) 



The equilibrium solution is obtained setting (p = Q. 

In silicon, where the peaks in the loss function are well 
defined, it is justified to use a single plasmon pole model 
V\;(t) = -iAfe{e-^"''^6'(r)-t-e*"''^6i(-T)} with plasmon 
energy ujk and intrinsic strengths A^ for each matrix ele- 
ment of W. Besides (p, the total exponent becomes then 
a-k [e*'^'='^ — l] with ak = \-kl^\ obtained from the cor- 
responding GW results. We find that a^ varies around 
0.3. Taylor expansion of the exponential leads then to 
the spectral function 



Akiuj) 



^ n! (w 

n=0 ^ 



efe + n(I)fc)^+r2/ 



(4) 



where efc — Re[efc] and Ffc = Im[efc]. Eq. ^ is sim- 
ilar to the plasmon pole version of the CE (cf. Ref. 
14|). However here the exponential solution arises from 
a straightforward approximation to the fundamental dif- 
ferential equation ([1]): the linearization of the Hartree 
potential reveals the boson of the model (i.e., the plas- 
mon via peaks in W), and the diagonal approximation of 
g gives rise to each isolated fermion. Our results are sum- 
marized in Fig. [TJ The dot-dashed line gives the result of 
this procedure together with the cross sections and the 
secondary electron background. The shapes of the QP 
peaks change little with respect to GW, but now the full 
series of satellites is present. The internal structure of 
the satellites which originate from the multiple valence 
bands, is also reproduced. This validates the decoupling 
approximation in the dense valence band region where, 
contrary to the case of an isolated core level, its suc- 
cess is a 'priori far from obvious. However, the intensity 
of the observed satellites is significantly underestimated. 



This discrepancy is similar to that found for the CE in 
simple metals, where extrinsic losses were suggested as 
a likely cause IJ]. These might also be reduced by in- 



terference effects [32|. To check this possibility we esti- 



mated the contributions from both effects to the satellite 
strengths a^ using Eq. (32) and (36) of Ref. [33[. This 
approach uses a plasmon pole model, Inglesficld fluctua- 
tion potentials, and an average over hole positions that 
takes account of the photoelectron mean free path A [33| . 
We observe that the averaged total satellite line shape in 
this model is similar to that for the intrinsic part, with 
a width 7 « 2 eV due to plasmon dispersion. Thus we 
can approximate the extrinsic and interference effects by 
renormalizing the intrinsic satellite intensity, i.e., by the 
replacement ak = Ofc + a^^* + «*"■'' in Eq. (jl)). These quan- 
tities are evaluated with ujp — 16.7 eV and A = 17.5 A at 
800 eV for Si, yielding a^^* = 0.63 and a"'-'' = -.11. This 
also modifies the strength Zk = e^"-'' of the QP peaks, 
but preserves overall normalization. The broadening of 
the satellites must also be increased, F -^ F -I- 7^7. The 
total spectrum thus obtained (black line) is in unprece- 
dented agreement with experiment. We stress that this 
result contains no fit parameters besides the two scaling 
factors (for spectrum and background) due to the arbi- 
trary units of the experiment. 

The success of our present approach stresses the need 
to go beyond the GWA. The exponential representa- 
tion of G implicitly corresponds to a vertex correction 
F = — ^'^- to the self-energy. Since our derivation yields 
Q as a function of the screened potential (p ([3]), this 
functional derivative can be performed explicitly, using 
.i£ll - r.-T^SGr.-i Prom Eq. ©, a straightforward 
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G-'§G-'- 



derivative of Q contains a series of satellite contributions. 
The two inverse Green's functions lead to a significant 
complication, because they contain the inverse of this se- 
ries. This clearly illustrates the difficulty of modelling 
F in order to treat dynamical effects. It suggests rather 
to concentrate on modelling j^, where the various con- 
tributions are simply summed, and hence to search for 
a self-energy of the form S = —iW^G^^ instead of 
S — iGWT. In conclusion, on the basis of our exper- 
imental XPS data we have analyzed the failure of GW 
to reproduce plasmon satellites and linked this failure 
to the appearance of a artificial plasmaron peak. On 
the other hand, GW results are fair when the imaginary 
part of S, hence the intensity of the corresponding plas- 
mon, is small enough so that no sharp plasmaron is cre- 
ated. Thus surprisingly, one might expect GW to work 
better in describing satellites stemming from local plas- 
mon or interband excitations close to the Fermi level in 
"strongly correlated" materials than for the strong plas- 
mon structures in conventional semiconductors. Starting 
from the fundamental equations of MBPT we have de- 
rived an exponential solution to the one-particle Green's 
function, analogous to that from the CE, that overcomes 



the drawbacks of the GWA. Comparison to new photoe- 
mission data shows that this yields a very good descrip- 
tion of the spectral function of bulk silicon, including the 
satellites series. By calculating the secondary electron 
background, cross section corrections as well as a cor- 
rection for extrinsic and interference effects, we achieve 
an agreement between theory and experiment that can 
be considered as a benchmark. Our derivation also sug- 
gests how the results can be improved in cases where the 
presently used approximations are inadequate. Finally, 
by accessing an expression for the vertex function, our 
approach yields precious hints for directions to take in 
modeling dynamical effects beyond the GWA. 
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